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We develop a theory of quadratic quantum measurements by a mesoscopic detector. It is shown 
that quadratic measurements should have non-trivial quantum information properties, providing, for 
instance, a simple way of entangling two non-interacting qubits. We also calculate output spectrum 
of a quantum detector with both linear and quadratic response continuously monitoring coherent 
oscillations in two qubits. 
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The problem of quantum measurements with meso- 
scopic solid-state detectors attracts considerable current 
interest (see, e.g., chapters on quantum measurements 
in [1]). This interest is motivated in part by important 
role of measurement in quantum computing, and in part 
by the possibility, provided by the mesoscopic structures, 
to study directly both in theory and experiment transi- 
tion between quantum and classical behavior in systems 
that are large on atomic scale. Although mesoscopic de- 
tectors can be quite different and include, e.g., quantum 
point contacts (QPC) [2-8], normal and superconduct- 
ing SET transistors [9-15], SQUID magnetometers [16] 
and generic mesoscopic conductors [17,18], the operat- 
ing principle of almost all of them is the same. Mea- 
sured quantum system controls transmission amplitude 
t of some particles (electrons, Cooper pairs, or magnetic 
flux quanta) between the two reservoirs, so that the flux 
of these particles provides information on the state of the 
system [19]. In a generic situation the amplitude t varies 
together with some control operator x, and for sufRciently 
weak detector-system coupling, the dependence t{x) can 
be approximated as linear. Dynamics of such linear mea- 
surements is well understood (see, e.g., [20,21]). 

At some special bias points, however, the linear re- 
sponse coefficient of the t{x) dependence vanishes and 
this dependence becomes quadratic. This can happen, 
for instance, if the amplitude t is formed by more than 
one interfering tunneling trajectories. Known examples 
of such situation include dc SQUIDs and superconduct- 
ing SET transistors. In this work, we show that quan- 
tum detector operating at such a special point should en- 
able measurements of product operators referring to sep- 
arate systems and have non-trivial quantum information 
processing properties, e.g., create simple entanglement 
mechanism for non-interacting qubits. Specifically, we 
consider the measurement of two qubits (Fig. 1) which is 
the simplest system that reveals non-trivial characteris- 
tics of quadratic detection. (Quadratic detectors can not 
measure individual qubits since dynamic variables of one 
qubit are given by the Pauli matrices a for which = 1.) 



The two qubits (indexed by j = 1,2) are assumed to be 
coupled to one detector through their basis-forming vari- 
ables, i.e., X = c\(j\ -I- c^cfi so that 
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(1) 



The last term in this equation appears due to non- 
linearity of t{x). If Sj = 0, A 7^ 0, one has purely 
quadratic detector. For the two qubits, Eq. (1) represents 
the most general dependence of t on , whereas for mea- 
surements of other systems, expansion in the measured 
operators similar to Eq. (1) can be justified as Taylor's 
expansion in weak detector-system coupling. 
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FIG. 1. Diagram of a mesoscopic detector measuring two 
qubits. The qubits modulate amplitude t of tunneling of de- 
tector particles between the two reservoirs. 



The Hamiltonian of the detector-qubit system is: 



(2) 



where Ho - -(1/2) ^,=1,2^^^ + A,aJ) + W2)a>2. 
Here Aj is the tunnel amplitude and Ej is the bias of the 
jth qubit, v is the qubit interaction energy, Hd is the de- 
tector Hamiltonian, and are the detector operators 
that create excitations when a particle is transferred, re- 
spectively, forward and backward between the detector 
reservoirs. For instance, for the QPC detector, ^t,^ de- 
scribe excitation of electron-hole pairs in the QPC elec- 
trodes. The qubit interaction is not essential to us, and 
we will frequently take v = below. 

We make two assumptions about the detector: the tun- 
neling between reservoirs is weak, so that the evolution of 
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the system can be described in the lowest non-vanishing 
order in t; characteristic time scale of tunneling is much 
shorter than that of the qubit evolution due to Hq- In 
the example of the QPC detector these assumptions mean 
that the QPC operates in the tunneling regime and that 
the voltage across it is much larger than the qubit ener- 
gies. Under these assumptions, precise form of the de- 
tector Hamiltonian is not important and dynamics of 
measurement is defined by the correlators 

7+= / dtm^}), 7-= / dt{£}m. (3) 
Jo Jo 

In Eq. (3), the angled brackets denote averaging over in- 
ternal degrees of freedom of the detector reservoirs which 
are taken to be in a stationary state. The correlators 
that do not conserve the number of par- 
ticles are assumed to vanish. 

The measurement contribution to the evolution of the 
qubit density matrix p is obtained by standard perturba- 
tion theory in tunneling and can be written down conve- 
niently in the "measurement" basis of eigenstates of the 
ai operators, | TT) , | Ti) > UT) ) and | J.J.). Each state \k) of 
this basis is characterized by the value tk of the transmis- 
sion amplitude (1): ti = io+X^j Sj+X, t2 = io+<5i— <52— A, 
H = to — 81+82 — A, t4 = to — J2j 8j+X. To describe qubit 
dynamics conditioned on particular outcome of measure- 
ment, we keep in the evolution equation the number of 
particles n transferred through the detector. Since the 
correlators that do not conserve n vanish, only terms di- 
agonal in n contribute to the evolution, and in the lowest 
order in tunneling, measurement contribution to p is: 

p^, = -(i/2)(r+ + r_)(|ifep + |t,|>^,+ 

r_tkttpl+' + T+t*M^' - i[8H, p^\u . (4) 

Here 8H. = fcjcr^ + 8v(t\(t\ is the renormalization 
of the qubit Hamiltonian due to coupling to the de- 
tector: fcj = Re((5jto + (5j'A*)Im(7_ -I- 7+) and 8v = 
Rei8i8*2 + toA*)Im(7_ + 7+), where j,f = 1,2, / ^ j. 
Equation (4) is the basis for our quantitative discussion 
of quadratic measurements. It generalizes to arbitrary 
detector and two qubits, the equation obtained in [5] for 
a qubit measured with the QPC in the tunnel regime. 

Disregarding the index n in Eq. (4), we obtain equa- 
tion for the measurement-induced evolution of the qubit 
density matrix averaged over different measurement out- 
comes. Together with the evolution due to the qubit 
Hamiltonian Hq this equation is: 

Pki = -IkiPki - i[Ho, p]ki ■ (5) 

Here jki = (l/2)(r+ + r_)|ife - t^p, with r± = 2Rej±, 
and we included in Hq two renormalization terms: 8H 

(4) and SH' due to phases (fki = arg(tfct;*) of the transfer 
amplitudes in Eq. (4) defined by: [8H',p]i.i = (r+ — 
r_)|tfcf(| sin (fikiPki- 



Evolution (5) of the qubit density matrix is reflected in 
the detector current. The form of the current / operator 
in the qubit space is obtained by the same lowest-order 
perturbation theory in tunneling that leads to Eq. (5): 

I={T+-T_)th. (6) 

This equation can be used to calculate both the dc cur- 
rent (/) = Trjlpo}, where po is the stationary solution 
of Eq. (5), and the current spectral density 

/•OO 

Si = So + 2 / drcosa;r(Tr{/e-^^[/po]} - {if). (7) 

Here So = {T+ + r_)TT{tHpo}, and e^''[A] denotes the 
evolution of the matrix A during time interval r governed 
by Eq. (5). 

Decay of the off-diagonal matrix elements Eq. (5) is 
the result of averaging over the measurement outcomes. 
However, since n is the classical detector output, it is le- 
gitimate to ask a question what is the qubit evolution for 
a specific measurement outcome n. Such a "conditional" 
description of the measurement dynamics (see, e.g., [21]) 
is convenient for calculation of more complicated correla- 
tors involved, for instance, in problems of feed-back con- 
trol of the dynamics of the measured system. For mea- 
surement process governed by Eq. (4) such a description 
is obtained by first solving this equation in terms of n. 
Noticing that Eq. (4) coincides in essence with the recur- 
rence relations for the modified Bessel functions /„ and 
assuming initial condition />fc;(0) = Pki{0)8nfi we get: 

Pliir) = pkimr+/T.r/Hn{2T\tkti\^T+rZ) ■ 

exp{-(l/2)(r+ + r_)(|ifc|2 + |ti|2)r - inw} ■ (8) 

The qubit density matrix conditioned on the particular 
"observed" number n of transferred particles is obtained 
then by selecting the term with this n in Eq. (8) and 
normalizing the resulting reduced density matrix. 

Quantitatively, conditional "Bayesian" equations for 
qubit evolution can be derived starting from Eq. (8) and 
following the same steps as in [21]. In particular, for weak 
detector-qubit coupling, |(5j|,|A| <C \to\, when individual 
tunneling events do not provide significant information 
on the qubit state, it is convenient to condition the evo- 
lution on the quasicontinuous current I(t) in the detector. 
Then, conditional equation for the evolution of the qubit 
density matrix is: 

Pki = -i[Ho,p]ki - jkipki 

+Ifit)pki Y^Pnih + I1- 2Ij) - i^ki) . (9) 

where So is the background current noise (see Eq. (7)), 
variation of which with the qubit state can be neglected 
in the weak-coupling limit, = (r_|_ -|-r_)|fo|^. Also, 
Ik = (r_|_ — r_) [ti; [^ is the average detector current in the 
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qubit state k, and If{t) = I(t) — pkklk is the fluctu- 
ation component of the detector current. Equation (9) 
is written in the Ito form, in which averaging over If{t) 
can be done by simply omitting the terms with it. In 
the weak-coupling regime, jki = (l/2)(r+ -|- r_)[(|t/c| - 
+ <^fe;|ioP]- It is the same ensemble-averaged de- 
cohcrcncc rate as in Eq. (5), but in Eq. (9), it leads to 
decoherence only after averaging over //(<). 

We now use equations obtained above to discuss sev- 
eral quantitative characteristics of quadratic measure- 
ments. We start with the purely quadratic detectors, 
when 5j = 0, so that 7i = 74 = (r+ - T_)\to + = /|| 
and I2 = h = (r+ — r„)|to — = J|j. In this case, 
if the qubits are stationary, Hq = 0, the detector effec- 
tively measures the product operator a^uf of the two 
qubits. I.e., on the time scale of measurement time 
Tm = 45'o/(/-ff — lu)"^^ the subspace {|1), |4)} in which 
the states of the two qubits are the same and the average 
detector current is is distinguished from the sub- 
space {|2), |3)} in which the states of the two qubits are 
opposite and the current is In, while the states within 
these subspaces are not distinguished. This property of 
quadratic measurements can be used to design simple 
error-correction scheme for dephasing errors [22] . 

Next, we consider the case of identical, unbiased, non- 
interacting qubits with non- vanishing Hamiltonian, Hq = 
— (A/2) a^. In this case the two degenerate zero- 
energy eigenstates of Hq can be chosen as {| ft) — | H 
)) I Ti) — I it)}- In the remaining subspace that will be 
denoted D+, in the basis {| TT) + | ii), | Ti) + I it)}, 
Hq reduces to —AcTx and mixes the states with similar 
and opposite states of the two qubits. Accordingly, there 
are three possible measurement outcomes characterized 
by the different dc currents (/) in the detector, 1^, 1^, 
and (/|| + I]-i)/2. These outcomes can be interpreted as 
measurement of the operator cr^cr^ + alal. Conditional 
equation (9) can be used to simulate how the qubits, on 
the time scale ~ ASq/I^, la = {1^ — I]i)/2, are driven 
into one of the three outcomes driven by the specific re- 
alization of the detector current. The probabilities of 
different outcomes depend on the initial qubit state. In 
the first two outcomes, the initial state is projected on 
one of the fully entangled states of the two qubits, e.g., 



{i)=in ^ IV') = (ln)-UT))/^/2. 



(10) 



Thus, quadratic measurements of two symmetric qubits 
provides a simple way of generating entangled states of 
qubits that in contrast to linear measurements [23] is 
based only on monitoring dc current instead of spectrum. 

In the third scenario, when (/) = [I^ + /||)/2, the 
two qubits are confined to the subspace D+ and perform 
coherent quantum oscillations. Equation for the density 
matrix (5) reduced to Dj^ is: 



Pkl = iA[(7a:,p]fei - 7 



0, P12 
P21 , 



where 7 = 2(r+ + r_)|Ap. Solving this equation and 
using the fact that the current operator (6) is reduced in 
D_^_ to /qCT^, we find the current spectral density (7): 



Si{uj) = 50-1- 



8A^7^7 



(a;2 _4A2)2 +^2^ 
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(12) 



where Sq = (r+ + r_)(|<oP + lAp). Qualitatively, for 
7 <C A, spectral density (12) describes coherent oscilla- 
tions of the two qubits with the frequency 2A, twice the 
oscillation frequency in one qubit. Similarly to the case 
of linear measurements [24] , the maximum of the ratio of 
the oscillation peak versus noise is 4. As one can see 
from Eq. (12), this maximum is reached in the case of 
weak measurement |A| <C |tol by the "ideal" detector for 
which arg(foA*) = and only r_|_ or r_ is non- vanishing. 

For different qubit tunnel amplitudes, transitions (with 
the rate (Ai — A2)^/27 for small Ai — A2) between the 
states I tt) ~ I ii) and | ti) — | it) niix the measurement 
outcomes 1^ and /| j . This means that for Ai 7^ A2 there 
are only two outcomes that have the same dc ciirrent and 
differ by current spectral densities. In one, the qubits 
are again in the subspace D+ and the spectral density 
is given by Eq. (12) where now 2A ^ Ai + A2. In the 
other, the qubit dynamics is confined to the subspace 
orthogonal to D+, the basis of which can be chosen as 
{|TT)-Ui), ITi)-|iT)}- Evolution of the qubit density 
matrix is described in this basis by the same Eq. (11) with 
A (Ai — A2)/2, and as a result the current spectral 
density is given by the same Eq. (12) with 2A ^ Ai — 
A2 and describes the current peak at difference of qubit 
oscillation frequencies. 

As the last application of the general theory we con- 
sider two identical qubits measured by a weakly and sym- 
metrically coupled detector with arbitrary non-linearity. 
It is convenient to discuss this situation in terms of the 
total effective spin S of the two qubits which determines 
the amplitude (1) of detector tunneling: 



t = tQ + + A(2S'2 - 1) , 



(13) 



The S = state (10) does not evolve in time un- 
der the qubit Hamiltonian and represents one of the 
measurement outcomes characterized by the dc detec- 
tor current 1^ and fiat current spectral density Si{u) = 
(r+ -|-r_)|to — Ap. Three other, 5=1, states are mixed 
by measurement and represent the second measurement 
outcome. We take the basis of the 5 = 1 subspace 
as three energy eigenstates 5^ = —1,0,1 with energies 
{A, 0, — A}. Interaction with the detector induces tran- 
sitions between these states with the rate independent 
of the transition's direction, so that the stationary qubit 
density matrix is po = 1/3. Equations (6) and (13) show 
then that the dc detector current for this outcome is: 



(11) (/) = (r+ - r_)[2(|io|^ + |A|2) + |io + A|2 + 8\5\']/3 , 
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and can be written as (/) = (/|| + In + /||)/3, where 
the currents /-fT'--- introduced in the same way as 
before, e.g, 7t| = (r+-r_)|to+25+Ap. The background 
current noise coincides with (7) with r+ — r_ replaced 

byr+ + r_. 

The system performs oscillations at frequencies A and 

2A whose spectral peaks should have Lorcntzian form 
for weak detector-qubit coupling. Evaluating the current 
matrix elements from Eqs. (6) and (13), and evolution 
of the density matrix from Eq. (5) reduced to the S = 1 
subspace, we obtain parameters of these Lorentzians in 
the spectral density (7): 



3(a;-iA)2+72 ' 



(14) 



J = 1,2, 7i = (r++r_)(|,5|2 + |A|V2), 72 = 271, 
ai = 4(r+ - r_)Re[(to + X)5*] = {1^ - hl)/'^ . 
a2 = 2(r+ - r_)(Re[toA*] + \5f) = {1^ + hi ' 2JTi)/4- 

Note that condition po = 1/3, used in Eq. (14) and Fig. 

2, implies that the coefficient S of linear measurement 
mixing all three 5=1 states, does not vanish identically. 
Otherwise, only two states are mixed and po = 1/2, as 
assumed in Eq. (12). In general, there is also a spectral 
peak at w = caused by switching between states with 
different average currents. For small but finite 5, 5 <C A, 
this peak can be very high, e.g., in the case of an "ideal" 
detector (r_ = argtoA* = argio^* = 0) its height and 
half-width are, respectively, {8X'^ /276'^)So and 6S'^r+. 

Figure 2 shows current spectral density calculated from 
Eqs. (5) and (7) for an ideal detector without making 
weak-coupling approximation. Figure 2a illustrates the 
transition between "single-qubit" oscillations at ~ A in 
the case of linear measurement and oscillations at w ~ 2 A 
for quadratic measurement. One can see that in agree- 
ment with Eq. (14) the ~ A peak is typically higher 
than the quadratic peak at double frequency. It is at 
the same time more sensitive to qubit-qubit interaction 
as illustrated in Fig. 2b. Even weak interaction <C A 
first suppresses and then splits the peak at w ~ A in two 
while affecting the quadratic peak only slightly. 

In summary, we discussed quadratic quantum measure- 
ments and have shown that they should have non-trivial 
properties providing, for instance, simple method of en- 
tangling non-interacting qubits. We also calculated out- 
put spectra of the quadratic detector measuring coherent 
oscillations in two qubits. Consistent with the case of 
classical oscillations, quadratic measurements results in 
the spectral peak at frequency that is twice the frequency 
of individual qubit oscillations. Quadratic measurements 
should be an interesting and potentially useful tool in 
solid-state quantum devices. 
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FIG. 2. Spectra of coherent quantum oscillation in two 
qubits measured with arbitrary non-linearity, (a) Transi- 
tion from linear to quadratic measurement in the case of 
non-interacting qubits, = 0, for r+|to|^ = 4. OA. Dashed 
line, solid line and dot-dashed line correspond to linear, 
{6 = O.lto, A = 0), intermediate {S = X = O.lto), and 
quadratic (5 = 0, A = O.lto) case, (b) Effect of qubit-qubit 
interaction on spectrum for r+ \to\'^ = 2. OA and S = \ = O.lto- 
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